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Construct a model of the received midamble 
sequence data, r, the channel response, ft , and 
the known midamble codes, such as 
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Reconfigure the model by replacing the known 
midamble codes with K right circulant matrix 
blocks, B, such as 



[M- 



Solve the reconfigure model using a least 
squares solution, such as by a discrete fourier 
transform or a single cyclic correlator solution 
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transmits K communication bursts in 
a shared spectrum in a time slot of a 
time division duplex communication 
system. The system associated with N 
midamble sequences. Each burst has an 
associated midamble sequence. A receiver 
receives a vector corresponding to the 
transmitted midamble sequences of the K 
communication bursts. A matrix having 
N identical right circulant matrix blocks is 
constructed based in part on the known N 
midamble sequences. The wireless channel 
between the transmitter and receiver is 
estimated based on in part one of the N 
blocks and the received vector. 
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CHANNEL ESTIMATION FOR TIME DIVISION 
DUPLEX COMMUNICATION SYSTEMS 

This application claims priority to U.S. Provisional Patent Application No. 
60/175,167, filed on January 7, 2000. 

5 BACKGROUND 

The invention generally relates to wireless communication systems. In 
particular, the invention relates to channel estimation in a wireless communication 
system. 

Figure 1 is an illustration of a wireless communication system 10. The 

10 communication system 10 has base stations 12, to 12 5 which communicate with user 
equipments (UEs) 14, to 14 3 . Each base station 12, has an associated operational 
area where it communicates with UEs 14, to 14 3 in its operational area. 

In some communication systems, such as code division multiple access 
(CDMA) and time division duplex using code division multiple access 

1 5 (TDD/CDMA), multiple communications are sent over the same frequency spectrum. 
These communications are typically differentiated by their chip code sequences. To 
more efficiently use the frequency spectrum, TDD/CDMA communication systems 
use repeating frames divided into time slots for communication. A communication 
sent in such a system will have one or multiple associated chip codes and time slots 

20 assigned to it based on the communication's bandwidth. 

Since multiple communications may be sent in the same frequency spectrum 
and at the same time, a receiver in such a system must distinguish between the 
multiple communications. One approach to detecting such signals is single user 
detection. In single user detection, a receiver detects only the communications from 

25 a desired transmitter using a code associated with the desired transmitter, and treats 
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signals of other transmitters as interference. Another approach is referred to as joint 
detection. In joint detection, multiple communications are detected simultaneously. 

To utilize these detection techniques, it is desirable to have an estimation of 
the wireless channel in which each communication travels. In a typical TDD system, 
5 the channel estimation is performed using midamble sequences in communication 
bursts. 

A typical communication burst 16 has a midamble 20, a guard period 18 and 
two data bursts 22, 24, as shown in Figure 2. The midamble 20 separates the two 
data bursts 22, 24 and the guard period 1 8 separates the communication bursts 1 6 to 

10 allow for the difference in arrival times of bursts 16 transmitted from different 
transmitters. The two data bursts 22, 24 contain the communication burst's data. The 
midamble 20 contains a training sequence for use in channel estimation. 

After a receiver receives a communication burst 16, it estimates the channel 
using the received midamble sequence. When a receiver receives multiple bursts 16 

15 in a time slot, it typically estimates the channel for each burst 16. One approach for 
such channel estimation for communication bursts 16 sent through multiple channels 
is a Steiner Channel Estimator. Steiner Channel Estimation is typically used for 
uplink communications from multiple UEs, 14, to 14^, where the channel estimator 
needs to estimate multiple channels. 

20 In some situations, multiple bursts 16 experience the same wireless channel. 

One case is a high data rate service, such as a 2 megabits per second (Mbps) service. 
In such a system, a transmitter may transmit multiple bursts in a single time slot. 
Steiner estimation can be applied in such a case by averaging the estimated channel 
responses from all the bursts 16. However, this approach has a high complexity. 

25 Accordingly, it is desirable to have alternate approaches to channel estimation. 

SUMMARY 

A single transmitter transmits K communication bursts in a shared spectrum 
in a time slot of a time division duplex communication system. The system 
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associated with N midamble sequences. Each burst has an associated midamble 
sequence. A receiver receives a vector corresponding to the transmitted midamble 
sequences of the K communication bursts. A matrix having N identical right 
circulant matrix blocks is constructed based in part on the known N midamble 
5 sequences. The wireless channel between the transmitter and receiver is estimated 
based on in part one of the N blocks and the received vector. 



BRIEF DESCRIPTION OF THE DRAWING(S) 
Figure 1 is a wireless communication system. 
Figure 2 is an illustration of a communication burst. 
10 Figure 3 is a simplified multiburst transmitter and receiver. 

Figure 4 is a flow chart of multiburst channel estimation. 



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT(S) 
Figure 3 illustrates a simplified multicode transmitter 26 and receiver 28 in a 
TDD/CDMA communication system. In a preferred application, such as a 2 Mbs 
15 downlink service, the receiver 28 is in a UE 14, and the transmitter 26 is in a base 
station 12,, although the receiver 28 and transmitter 26 may be used in other 
applications. 

The transmitter 26 sends data over a wireless radio channel 30. The data is 
sent in K communication bursts. Data generators 32, to 32 K in the transmitter 26 

20 generate data to be communicated to the receiver 28. Modulation/spreading and 
training sequence insertion devices 34, to 34 K spread the data and make the spread 
reference data time-multiplexed with a midamble training sequence in the appropriate 
assigned time slot and codes for spreading the data, producing the K communication 
bursts. Typical values of K for a base station 12, transmitting downlink bursts are 

25 from 1 to 16. The communication bursts are combined by a combiner 48 and 
modulated by a modulator 36 to radio frequency (RF). An antenna 38 radiates the 
RF signal through the wireless radio channel 30 to an antenna 40 of the receiver 28. 
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The type of modulation used for the transmitted communication can be any of those 
known to those skilled in the art, such as binary phase shift keying (BPSK) or 
quadrature phase shift keying (QPSK). 

The antenna 40 of the receiver 28 receives various radio frequency signals. 
5 The received signals are demodulated by a demodulator 42 to produce a baseband 
signal. The baseband signal is processed, such as by a channel estimation device 44 
and a data detection device 46, in the time slot and with the appropriate codes 
assigned to the transmitted communication bursts. The data detection device 46 may 
be a multiuser detector or a single user detector. The channel estimation device 44 

10 uses the midamble training sequence component in the baseband signal to provide 
channel information, such as channel impulse responses. The channel information 
is used by the data detection device 46 to estimate the transmitted data of the received 
communication bursts as hard symbols. 

To illustrate one implementation of multiburst channel estimation, the 

15 following midamble type is used, although multiburst channel estimation is 

applicable to other midamble types. The K midamble codes, m {k) , where k = 1 ,..K, 
are derived as time shifted versions of a periodic single basic midamble code, m p , 
of period P chips. The length of each midamble code is L m = P + W- 1 . W is the 
length of the user channel impulse response. Typical values for L m are 256 and 5 1 2 

20 chips. W is the length of the user channel impulse response. Although the following 
discussion is based on each burst having a different midamble code, some midambles 
may have the same code. As, a result, the analysis is based on N midamble codes, 
N < K. Additionally, the system may have a maximum number of acceptable 
midamble codes N. The receiver 28 in such a system may estimate the channel for 

25 the N maximum number of codes, even if less than N codes are transmitted. 

The elements of m p take values from the integer set { 1 , - 1 } . The sequence m )t 

is first converted to a complex sequence m P [i] = f *m P [i], where / = V -P . The 
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m m are obtained by picking K sub-sequences of length L m from a 2^ long 

sequence formed by concatenating two periods of Wp . The /'* element of m {k) is 

related to m p by Equation 1 . 

m\ k) = m P [(K - k)W + /], for 1 < i < P - (K - k)W 

= m P [i-P+(K-k)w] f foxP-{K-k)W<i<P + W-\ 

Equation 1 

Thus, the starting point of m {k \ k = l..K shifts to the right by W chips as k increases 
from 1 to K. 

The combined received midamble sequences are a superposition of the K 
convolutions. The K h convolution represents the convolution of m (k) with . 

Jft) is the channel response of the kf h user. The preceding data field in the burst 

corrupts the first (W - 1) chips of the received midamble. Hence, for the putpose of 

channel estimation, only the last P of L m chips are used to estimate the channel. 

Multiburst channel estimation will be explained in conjunction with the flow 
chart of Figure 4, To solve for the individual channel responses , Equation 2 is 
used. 
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r w ~-r LM are the received combined chips of the rnidamble sequences. The m values 
are the elements of m p . 

Equation 2 may also be rewritten in shorthand as Equation 3. 



Equation 3 



k-\ 



Each M {k) is a KW-by-W matrix, r is the received rnidamble chip responses. 

When all the bursts travel through the same channel, h {]) ->-h {k) can be replaced by h 
as in Equation 4, 50. 



h=f 



Equation 4 



G is defined as per Equation 5. 

G e [M°\..., A/ ( *>,..., M lK) j Equation 5 

As a result, G is a KW-by-KW matrix. Since G is a right circulant matrix, Equation 
4 can be rewritten using K identical right circulant matrix blocks B, as per Equation 
6, 52. 

' B 



!><*> 



= D 



Equation 6 



B is a W-by-W right circulant matrix. The number of B-blocks is K. Using Equation 
6, Equation 4 can be rewritten as Equation 7. 



Dh =r 



Equation 7 
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Equation 7 describes an over-determined system with dimensions KW-by-W, One 
approach to solve Equation 7 is a least squares solution, 54. The least squares 
solution of Equation 7 is given by Equation 8. 



h = (d h d)~ } D H F 



Equation 8 



D H is the hermitian of D. 

Applying Equation 6 to Equation 8 results in Equation 9. 



{d»dY = \{b»b)- 



Equation 9 



The received vector r of dimension KW can be decomposed as per Equation 10. 



Equation 10 



10 The dimension of r k is W. Substituting Equations 9 and 10 into Equation 8, the 
least-squares solution for the channel coefficients per Equation 1 1 results. 



Equation 1 1 



r k represents the average of the segments of r . Since Bis a square matrix, Equation 
11 becomes Equation 12. 



15 



h = B~*n 



Equation 12 
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Since B is a right circulant matrix and the inverse of a right circulant matrix is also 
right circulant, the channel estimator can be implemented by a single cyclic correlator 
of dimension 57, or by a discrete fourier transform (DFT) solution. 

A W point DFT method is as follows. Since B is right circulant, Equation 1 3 
can be used. 



Dp is the W point DFT matrix as per Equation 14. 



Equation 13 



10 
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Equation 14 



A c is a diagonal matrix whose main diagonal is the DFT of the first column of B, 
as per Equation 15. 



A c =diag(D w (B^\))) 



Equation 15 



15 



W = e' jW . Thus, D w is the DFT operator so that D w x represents the W point DFT 

of the vector x . By substituting Equation 13 into Equation 12 and using dz) - ^ 

« w 



results in Equation 16. 



Equation 16 
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D* w is the element-by-element complex conjugate of D w . 

Alternately, an equivalent form that expresses ft in terms of A R instead of 
A c can be derived, is a diagonal matrix whose main diagonal is the DFT of the 
first row of B per Equation 17. 

= diag(p w {B(\,tyj Equation 17 

Since the transpose of B, B T , is also right circulant and that its first column is the first 
row of B, B T can be expressed by Equation 1 8. 

B T - D~y) • h H -D w Equation 1 8 

Using Equation 18 and that Dj v = D w , A r R = A^ and that for any invertible matrix 

A, ^A T } 1 - {^) > B can be expressed as per Equation 19. 

B = D w • A R ■ D^} Equation 1 9 

Substituting Equation 19 into Equation 12 and that D^ 1 = results in Equation 20. 

W 




Equation 20 



Equations 16 or 20 can be used to solve for h . Since all DFTs are of length W, the 
complexity in solving the equations is dramatically reduced. 

An approach using a single cycle correlator is as follows. Since B~ ] is the 
inverse of a right circulant matrix, it can be written as Equation 21 . 
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Equation 21 
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The first row of the matrix T is equal to the inverse DFT of the main diagonal of A"^ . 
Thus, the matrix T is completely determined by A"^ 1 . 

The taps of the channel response h are obtained successively by an inner 
product of successive rows of T with the average of W-length segments of the 
received vector F . The successive rows of T are circularly right shifted versions of 
the previous row. Using registers to generate the inner product, the first register 
holds the averaged segments of F , and the second register is a shift register that holds 
the first row of the matrix T. The second register is circularly shifted at a certain 
clock rate. At each cycle of the clock, a new element of h is determined by the 
inner product of the vectors stored in the two registers. It is advantageous to shift the 
first row of the matrix T rather than the received midambles. As a result, no extra 
storage is required for the midambles. The midambles continue to reside in the 
received buffer that holds the entire burst. Since the correlator length is only W, a 
significant reduction in complexity of estimating the channel is achieved. 
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CLAIMS 

What is claimed is: 

L A method for estimating a wireless channel in a time division duplex 
communication system using code division multiple access, the system associated 
with N midamble sequences, the wireless channel existing between a single 
transmitter and a single receiver, the single transmitter transmitting K communication 
5 bursts in a shared spectrum in a time slot, each burst having an associated midamble 
sequence of the N sequences, the receiving knowing the N midamble sequences, the 
method comprising: 

receiving a vector corresponding to the transmitted midamble sequences of the 
K communication bursts at the single receiver; 
1 0 constructing a matrix having N identical right circulant matrix blocks based 

in part on the known N midamble sequences; and 

estimating the wireless channel based on in part one of the N blocks and the 
received vector. 

2. The method of claim 1 wherein the wireless channel estimating is 
performed using a least squares solution, 

3. The method of claim 2 wherein the least squares solution is 
implemented using a single cyclic correlator. 

4. The method of claim 2 wherein the least squares solution is 
implemented using a discrete fourier transform solution. 

5. The method of claim 1 wherein N is a maximum number of midamble 
codes associated with the system. 
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6. The method of claim 1 wherein N is a number of different midambles 
transmitted in the K bursts. 

7. The method of claim 1 wherein a dimension of the estimation is a 
length of a channel response of the wireless channel. 

8. A receiver for use in a wireless time division duplex communication 
system using code division multiple access, the system associated with N midamble 
sequences, a single transmitter in the system transmits K communication bursts in a 
shared spectrum in a time slot, each burst having an associated midamble sequence 

5 of the N sequences, the receiver knowing the N midamble sequences, the receiver 
comprising: 

an antenna for receiving the K communication bursts including a vector 
corresponding to the transmitted midamble sequences of the bursts; 

a channel estimator for constructing a matrix having N identical right 
10 circulant-matrix blocks based in part on the known N midamble sequences and 
estimating the wireless channel between the receiver and the single transmitter based 
on in part one of the N blocks and the received vector; and 

a data detector for recovery data from the received communication bursts 
using the estimated wireless channel. 

9. The receiver of claim 8 wherein the data detector is a multiuser 
detector. 

10. The receiver of claim 8 wherein the data detector is a single user 
detector. 

11. The receiver of claim 8 wherein the wireless channel estimating is 
performed using a least squares solution. 
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12. The receiver of claim 11 wherein the least squares solution is 
implemented using a discrete fourier transform solution. 

13. The receiver of claim 11 wherein the least squares solution is 
implemented using a single cyclic correlator. 

14. A wireless spread spectrum communication system using code division 
multiple access associated with N midamble sequences, the system communicating 
using communication bursts, each burst having an associated midamble sequence, the 
system comprising: 

5 a base station comprising: 

a data generator for generating data; 

a plurality of modulation/spreading devices for formatting the generated 
data into K communication bursts time multiplexed to be in a same time slot and in 
a shared spectrum; and 
10 an antenna for radiating the K communication bursts; and 

a user equipment comprising: 

an antenna for receiving the K communication bursts including a vector 
corresponding to the transmitted midamble sequences of the bursts; 

a channel estimator for constructing a matrix having N identical right 
1 5 circulant matrix blocks based in part on the N midamble sequences and estimating 
the wireless channel between the base station and the user equipment based on in part 
the K block matrix and the received vector; and 

a data detector for recovering data from the received communication bursts 
using the estimated wireless channel, 

15. The system of claim 14 wherein the data detector is a multiuser 
detector. 
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1 6 . The system of claim 1 4 wherein the data detector is a plurality of single 
user detectors. 

17. The system of claim 14 wherein the wireless channel estimating is 
performed using a least squares solution. 

18. The system of claim 17 wherein the least squares solution is 
implemented using a discrete fourier transform solution. 

19. The system of claim 17 wherein the least squares solution is 
implemented using a single cyclic correlator. 

20. The system of claim 14 wherein the base station is effectively 
transmitting data at a 2 Mbps data rate to the user equipment. 
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